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The shoulder structure of charged particles multiplicity distributions (MD's) 
^ ! in full phase space in e"'"e~ annihilation at the Zq peak and the quasi-oscillatory 

' behavior of the ratio of factorial cumulants over factorial moments, ifg, as a 

function of the order g, are quantitatively reproduced within a simple parametri- 
zation of the MD in terms of a weighted superposition of two Negative Bino- 
mial Distributions, associated to two- and multi-jet production, i.e., to hard 
gluons radiation. 
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1. Introduction 



A definite understanding of tlie multiparticle production process in full phase space in 
e^e" annihilation at the Zq peak is still lacking. Experimental data show two interesting 
features: a clear shoulder is visible in the intermediate multiplicity range in charged 
particles multiplicity distributions (MD's)|l], ^, ^; the ratio of factorial cumulants over 
factorial moments of the MD, Hq, when plotted as a function of its order g, decreases 
sharply to a negative minimum at g=5 and follows then a quasi-oscillatory behavior ^. 
In order to interpret these results, two complementary approaches can be followed. First, 
properties of multiparton final states can be computed in the framework of perturbative 
QCD by exploiting its branching structure[P, 0]; partonic and hadronic distributions are 
then directly compared by assuming the Local Parton Hadron Dualityj^ or its generalized 
version as hadronization prescription. Second, a phenomenological approach can be 
taken and data compared with simple parametrizations, like for instance the Negative 
Binomial (NBD)|lOl or the log-normal (LND)[0 distributions. 

Within the first approach, the partons' MD has been computed in Double Log Ap- 
proximation (DLA)I^, getting a much wider MD than the experimental one. By adding 
an approximate treatment of the energy-momentum conservation law||12||, the predicted 
MD becomes narrower at non asymptotic energies, but it can quantitatively reproduce 
neither the general shape of the experimental MD nor its shoulder. Concerning the ratio 
Hq, QCD predictions beyond DLA have been considered [0]. They qualitatively show the 
same behavior of experimental data; the result supports the general picture of LPHD, but 
it is quite far from being quantitatively satisfactory. 

Within the second more phenomenological approach, let us consider the two most 
popular parametrizations of MD's, the NBD and the LND. Both of them are related to 
QCD. The NBD is the exact solution for the gluon MD in a quark-jet in Leading Log 



Approximation |T^ and more detailed QCD calculations in DLA plus recoil corrections 
turn out to be close to NB behavior |T^. It reveals the self-similar and Markoffian 
structure of parton shower evolution; it is usually interpreted in the framework of the clan 



model |]10[ as an indication of the two-step nature of the production process. The LND 



results from the application of the central-limit theorem to self-similar parton showers|ll 



Both parametrizations quantitatively describe e^e annihilation data in full phase space 
at cm. energies below the Zq peak[|l^, |1^ and qualitatively reproduce the overall shape 
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of the MD at the Zq peak, but neither the NBD nor the LND reproduces the shoulder 
structure, although the LND fares generally better in terms of x^/NDF. Concerning the 
ratio Hq of factorial cumulants over factorial moments, parametric models give a poor 
description of experimental results. A full NBD predicts indeed positive definite and 
monotonically decreasing Hq's. If one truncates a NBD in the high multiplicity tail to 
properly take into account the experimental cut due to finite statistics, one recovers a 



quasi-oscillatory behavior qualitatively consistent with experimental results[0; the size of 
the experimental effect is however strongly underestimated 0. Predictions of a truncated 
LND are in better agreement with data[Q, even though they do not provide a quantitative 
description of the experimental effect. 

In this letter, we propose a simple phenomenological parametrization of the MD in 
e~^e^ annihilation in full phase space in terms of a weighted superposition of two NBD's; 
both the shoulder structure of the MD and the quasi-oscillatory behavior of the ratio Hg 
are quantitatively reproduced. These results suggest to explain both effects in terms of 
the superposition of samples of events with a fixed number of jets, where the standard 
NB parametrization is indeed recovered. 



2. The shoulder structure 



The DELPHI Collaboration has shown in ITOll that the shoulder structure in the MD in 



e'^e^ annihilation can be explained by the superposition of the MD's coming from events 
with 2, 3 and 4 jets, as identified by a suitable jet-finding algorithm, and that the MD's in 
these classes of events are well reproduced by a single NBD for every value of the jet-finder 
parameter, i/cut- It is thus suggested that the shoulder is associated with the radiation 
of hard gluons resulting in the appearance of one of more extra jets in the hadronic final 
state. 

Let us also remind that a shoulder structure similar to the one observed in e+e~ 



annihilation has been observed in pp collisions at high energies ||2CI|| and was shown to be 
well described by a 5-parameter parametrization in terms of the weighted superposition 
of two NBD'sl21[. 



Inspired by these experimental results and NB universality for all classes of reactions, 
we propose now to parametrize experimental data on MD's in full phase space in e+e~ 
annihilation at the Zq peak in terms of the superposition of 2 NBD's; according to DEL- 
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phi's result, we argue that the two MD's should be associated with the contribution of 
two- and multi(> 3)-jets events respectively. We fix therefore the weight parameter to be 
equal to the fraction of 2-jet events. Formally, we perform a fit on available experimental 
data obtained at LEP|1|, ^ and SLC[^ with the following 4-parameter distribution: 

Pnia; ni, h; ^2, ^2) = «P„^^(ni, h) + (1 - a)Pi'^(n2, (1) 

where P^^{n,k) is the NBD, expressed in terms of two parameters, the average multi- 
plicity n and the parameter k, linked to the dispersion by D'^/n'^ = l/fi + 1/k, as: 

r>NBt- ,\ k{k + 1) . . . {k + n - 1) ( k \^ f n 



n,k)= ^ ^ --T --T (2) 

nl \n + k J \n + k / 

The weigth a gives the 2-jet events fraction: since the latter depends on the jet-finder 
parameter, i/cut, different values of a corresponding to different values of i/cut have been 
considered. 

As far as MD's in full phase space are concerned, one has also to take care of the 
"even-odd" effect, i.e., of the fact that the total numer of final charged particles must be 
even due to charge conservation; accordingly, the actual form used in the fit procedure is 
given by: 

_ j APn{a]ni,ki;n2,k2) if n even 

1 otherwise 
where A is the normalization parameter, so that J2'^=o-Pn = 1- 

The parameters of the fits and the corresponding x^/NDF for different values of a 

are given in Table 1; for comparison, the parameters of the NBD's fitted to 2- and 3-jet 

data samples by DELPHI Coll. |[T^ are also shown. As an example. Figure 1 compares the 

experimental MD's with the parametrization (|^) with best-fit parameters given in Table 1 

for a = 0.767; the two NBD contributions are plotted with a dotted line. The residuals, 

i.e., the difference between data and theoretical predictions, are also shown below each 

MD in units of standard deviations. One concludes that the proposed parametrization 

can reproduce the experimental data very well; no structure is visible in the residuals. 

The agreement holds for all considered values of a, even though the parameters of the 

two NBD's obviously depend on the choice of the weight parameter. To further support 

the interpretation of the two components of eq. (|l]) as the two- and multi-jet contribution 

respectively, we notice that the parameters of the two NBD's extracted from the fit at 

^Data on MD by L3 Coll.|22| have not been analyzed due to their large systematic error. 



4 



a given a are close to the fitted values extracted by DELPHI Coll. for 2- and 3-jet 
MD's at the corresponding value of ycuu in fact for these values of ycut the contribution 
of (> 4)-jet events is negligible (less than 10%). Finally, let us also notice that the fit 
parameters present large errors; their determination would be strongly improved by direct 
experimental analyses which could take properly into account the full covariance matrix 
of the MD's. These analyses are eagerly awaited. 

3. The ratio Hq 

Let us consider now the ratio of factorial cumulant over factorial moments 

a, = KJF, (4) 

as a function of the order q. The factorial moments, Fg, and factorial cumulant moments, 
Kq, can be obtained in general from the MD, P„, through the relations: 

oo 

Fq = Y,n{n-l)...{n-q + (5) 

n=q 

and 

Kq = Fq-Y.^(^'^^^Kq_,F,. (6) 

A direct experimental analysis of the ratio with correlations and systematics effects 
fully taken into account has been so far performed by SLD Coll.0] only. To extend 
the analysis to LEP experiments, we extracted a posteriori the values of the Hq from 
published MD's. Due to the lacking of the full covariance matrix for MD's, the statistical 
error on Hq cannot be properly calculated; in order to estimate it, for each experimental 
MD we calculated the ratio Hq from 1000 MD's obtained from the original one by allowing 
gaussian fluctuations around the measured value of each Pn the width of the gaussian being 
given by the experimental statistical error. The standard deviation of the Hq distribution 
for each order q was then taken as an estimate of the error associated with the ratio Hq 
itself. As a cross-check, this method has been applied to SLD data, since both MD and 
the ratio Hq are available in this case; a good agreement between the errors extracted 
with this procedure and the sum of statistical and systematic errors measured by SLD has 
been found. It should be stressed of course that this method can give just an approximate 
value of the size of the errors and direct experimental analyses by LEP Collaborations are 
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awaited. Let us now look at the ratio Hg as predicted by the parametrization proposed 
in Section 2. Since the Hg^s were shown to be sensitive to the truncation of the tail due 



to the finite statistics of data samples |T8[], moments were extracted from a truncated MD 
defined as follows: 

ptrunc _ f ^'Pnia; Ui, ki, n2, k2) if {rimin <n< Umax) and n even 
1 otherwise 

Here n^m and Umax are the minimum and the maximum observed multiplicity and A' is 
a new normalization parameter^, so that I]n=n'li„ p^™^<^ = i. 

Figure 2 shows the ratio Hg as a function of the order q for SLD]^, ALEPH|0, 
DELPHIfll] and OPAL[gl. It should be noticed that SLD, DELPHI and OPAL give very 
similar results, whereas ALEPH shows a different behavior, in particular at low orders 
of q and after, when oscillations start. However, due to the large errors and to the 
uncertainties in their estimate, a definite answer on the consistency among different set of 
data requires a direct experimental determination of errors. In Figure 2 predictions of the 
parametrization (|^ with parameters fitted to reproduce the MD as given in Table 1 for 
different values of a are also shown (solid lines). It turns out that this parametrization 
reproduces not only the shoulder observed in experimental MD's but also quantitatively 
describes the experimental behavior of the ratio Hg. Results are essentially independent 
of a; the small spread among theoretical predictions for different values of a can be taken 
as an estimate of the theoretical error. The plot referring to SLD shows in addition 
predictions of eq. @, i.e., of the same parametrization, but without taking into account 
the effect of truncation (dashed lines): the quasi-oscillations become smaller. To further 
support the interpretation of this oscillatory behavior as due to the superposition of two- 
and multi-jet events, it should be noticed that oscillations in samples of isolated 2- and 
3-jet events are much smaller than in the full sample of events[^]. In conclusion, the 
observed behavior of i/g's results from the convolution of two different effects, a physical 
one, i.e., the superposition of two components, and a statistical one, i.e., the truncation 
of the tail due to the finite statistics of data samples. 



^We found that a small (0.1%) error in the normalization leads to very different values for the i/^'s 
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4. Conclusions 

Two experimental features of multiparticle production in full phase space in e^e~ annihi- 
lation at the Zq peak, i.e., the shoulder in the MD and the oscillatory behavior of the ratio 
of factorial cumulants over factorial moments, Hq, as a function of the order have been 
adressed: a simple phenomenological parametrization of the MD in terms of a weighted 
superposition of two components, each one being distributed according to a NBD, has 
been proposed. The weight of the first component was taken equal to the fraction of 2-jet 
events, i.e., the two components were identified with the two- and the multi-jet contri- 
bution respectively. The shoulder structure is found to be quantitatively reproduced by 
this parametrization; the behavior of the ratio H,^ is also quantitatively described, after 
taking properly into account the effect due to the truncation of the tail. A consistent 
picture seems then to emerge: both experimental effects which deviate by a single NBD 
behavior are explained in terms of a common mechanism, linked to the emission of hard 
gluons. The simple NB parametrization is reestablished at the level of events with fixed 
number of jets. Further tests of the above mentioned picture can be provided by direct 
experimental analyses which take into account the effects of correlations and systematics; 
these analyses are eagerly awaited. 
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6. Table Captions 

Table 1. Parameters and xV^DF of the fit to experimental data from ALEPH0, DELPHIll|], 
OPAL[0 and SLD0] with the weighted sum of two NBD's. Results are shown for 
different values of a corresponding to the fraction of 2-jet events, /, experimentally 
measured by DELPHI Collaboration at different values of the jet-finder parameter 
ycut- NBD parameters extracted by the DELPHI Collaboration by fitting MD's of 
samples of events with 2- and 3-jets at different values of ycut are also shown for 
comparison in the last columns. 



7. Figure Captions 

Figure 1. Charged particles' MD in full phase space, P„, at the peak from ALEPH|p, 
DELPHI [|, SLD|] and OPAL|] are compared with equation (D with a = 0.767 
(see Table 1 for the values of the corresponding parameters) (solid lines); dotted lines 
indicate the two separate NBD contributions. The lower part of the figure shows the 
residuals, i.e., the difference between data and theoretical predictions, expressed 
in units of standard deviations. 

Figure 2. The ratio of factorial cumulnats over factorial moments, Hq as a function of q\ 
experimental data (diamonds) are compared with equation (^ for different values 
of a as in Table 1 (solid lines). SLD, ALEPH, DELPHI and OPAL data are explicitly 
indicated. The dashed lines in the first plot (SLD) show predictions of equation (^, 
i.e., the same parametrization as the solid lines but without taking into account the 
effect of truncation. In the figure only statistical errors of SLD data are shown. 
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Table 1 





ALEPHlD 


DELPHI [0 


OPALg] 


SLDi 


DELPHI 119 




a = 0.463 


y,ut = 0.02 /=0.463 


ni 

h 

n2 
k2 

XVNDF 


17.7±1.1 
111±168 
23.6±0.8 
32±15 
3.56/22 


18.2±0.2 
90±20 

23.9±0.2 

31±3 
8.95/21 


18.4±0.2 
71±11 

24.0±0.2 

28±2 
3.32/21 


18.4±0.2 
47±4 

23.0±0.2 
29±2 
17.6/21 


n2-jet 
k2-jet 

k^-jet 


18.5±0.1 

57±4 
22.9±0.1 
44±2 




a = 0.659 


y,ut = 0.04 /=0.659 


ni 
ki 
n2 
k2 

XVNDF 


18.5±0.7 
66±46 

25.5±1.0 
47±33 

3.72/22 


18.9±0.2 
63±8 

25.8±0.3 
44±5 

10.1/21 


19.0±0.1 
54±5 

25.9±0.2 
40±5 

4.40/21 


18.9±0.1 
42±3 

24.7±0.2 
37±3 

16.3/21 


n2-jet 
ks-jet 


19.4±0.1 
44±2 

24.8±0.1 
42±2 




a = 0.767 


Vcut = 0.06 /=0.767 


ni 
ki 
n2 
k2 

XVNDF 


19.1±0.5 
53±24 

27.0±1.1 
65±62 
3.86/22 


19.4±0.2 
52±6 

27.3±0.3 
61±10 
11.7/21 


19.5±0.07 

46±3 
27.5±0.2 

55±8 
6.30/21 


19.3±0.09 
39±2 

26.0±0.2 
47±5 
15.6/21 


n2-jet 
k2-jet 

k^-jet 


20.0±0.1 

38±1 
26.0±0.1 

45±2 




a = 0.834 


Vcut = 0.08 /=0.834 


ni 
ki 
n2 
k2 

XVNDF 


19.5±0.4 

45±15 
28.2±1.2 
92±121 
3.99/22 


19.8±0.1 
46±3 

28.6±0.3 
85±18 
13.9/21 


19.9±0.6 

40±2 
28.8±0.2 

76±15 
8.81/21 


19.6±0.1 
37±2 

27.1±0.3 
59±7 

15.2/21 


n2-jet 
k2-jet 

ks-jet 


20.4±0.1 
34±1 

26.8±0.1 
49±1 
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